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Abstract 

The behaviour of the magnetic field of a neutron star with a supercon- 
ducting quark matter core is investigated in the framework of the Ginzburg- 
Landau theory. We take into account the simultaneous coupling of the diquark 
condensate field to the usual magnetic and to the gluomagnetic gauge fields. 
We solve the Ginzburg-Landau equations by properly taking into account the 
boundary conditions, in particular, the gluon confinment condition. We found 
the distribution of the magnetic field in both the quark and hadronic phases 
of the neutron star and show that the magnetic field penetrates into the quark 
core in the form of quark vortices due to the presence of Meissner currents. 
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1 Introduction 



Recently, possible formation of diquark condensates in QCD at finite density has 
been re- investigated in series of papers following Refs.[0],0. It has been shown 
that in chiral quark models with non-perturbative 4-point interaction motivated 
from instantons or non-perturbative gluon propagators 0], P the anomalous 
quark pair amplitudes in the color antitriplet channel can be very large: of the 
order of 100 MeV. Therefore, one expects the diquark condensate to dominate the 
physics at densities beyond the deconfinement /chiral restoration transition density 
and below the critical temperature (of the order of 50 MeV). Various phases are 
possible. The so called two-flavor (2SC) and three-flavor (3SC) phases allow for 
unpaired quarks of one color. It has been also found f^, that there can exist a 
color-flavor locked (CFL) phase for not too large strange quark masses |p, where 
color superconductivity is complete in the sense that diquark condensation results 
in a pairing gap for the quarks of all three flavors and colors. 

The high-density phases of QCD at low temperatures are relevant for the ex- 
planation of phenomena in rotating massive compact stars which might manifest 
themselves as pulsars. Physical properties of these objects (once being measured) 
could constrain our hypotheses about the state of matter at the extremes of den- 
sities. In contrast to the situation for the cooling behaviour of compact stars 0, 
[ p!0|] , where the CFL phase is dramatically different from the 2SC and 3SC phases. 



we don't expect qualitative changes of the magnetic field structure for these phases. 
Therefore, below we shall restrict ourselves to the discussion of the simpler two- 
flavor theory first. Bailin and Love []11| used a perturbative gluon propagator which 
yielded a very small pairing gap and they concluded that quark matter is a type I 
superconductor, which expells the magnetic field of a neutron star within time-scales 
of 10^ years. If their arguments would hold in general, the observation of life-times 
for magnetic fields as large as lO^years |T2[, []T3| would exclude the occurence of 



an extended superconducting quark matter core in pulsars. These estimates are not 
valid for the case of diquark condensates characterized by large quark gaps. Besides, 
in Ref. |T^ the authors found that within recent non-perturbative approaches for 



the effective quark interaction that allow for large pairing gaps the quark condensate 
forms a type II superconductor. Consequently for the magnetic field H < Hd there 
exists a Meissner effect and for Hc2 > H > Hd the magnetic field can penetrate 
into the quark core in quantized flux tubes. However, they have not considered 
in that paper the simultaneous coupling of the quark fields to the magnetic and 
gluomagnetic gauge fields. 

Though color and ordinary electromagnetism are broken in a color supercon- 
ductor, there is a linear combination of the photon and the gluon that remains 
massless. The authors of Ref. [0 have considered the problem of the presence of 
magnetic fields inside color superconducting quark matter taking into account the 
possibility of the so called "rotated electromagnetism". They came to the conclusion 
that there is no Meissner effect and the external static homogeneous magnetic field 
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can penetrate into superconducting quark matter because in their case it obeys the 
sourceless Maxwell equations. To our opinion this result is obtained when one does 
not pose correct boundary conditions for the fields. Obviously it is energetically 
favorable to expell the magnetic field rather than to allow its penetration inside the 
superconducting matter. Using for the description of the diquark condensate inter- 
acting with two gauge fields the same model as in Refs. p|, |], |T5|, , the authors of 
Ref. [|17| have shown that the presence of the massless excitation in the spectrum 
does not prevent the Meissner currents to effectively screen static, homogeneous, 
external magnetic fields. In the present paper we will extend those studies to the 
consideration of inhomogeneous, vortex-type external fields. 



In Ref. |TH] we have derived the Ginzburg-Landau equations of motion for the 



diquark condensate placed in static magnetic and gluomagnetic fields 

e ^ q 
3 V3 



adp + p{dpd;)dp + 7(i V A +^ G8)^dp = 0, (1) 



where dp is the order parameter, a = t dn/dE, (3 = (dn/dE) 7 C(3)(7rTc)^^/8, •y = Pp 
/3/(6/i^), dn/dE = pp /i/vr^, t = (T — Tc)/Tc, Tc being the critical temperature, pp- 
the quark Fermi momentum, and for the gauge fields 



2 2 ^ sma{dpVdl- dlVdp) . ^ 

A„rot rot A + sm a A= i ^ — h sm a cos a Gs (2) 

^ 2qd^ 



2 ^ 2 ^ cos a{dpS7 dl — dlV dp J 

A.rot rot Gs + cos a G8= —i ^ — 1- sin a cos a A ■ (3) 

^ 2qd^ 



These equations introduce a "new" charge of the diquark pair q = A/r/^e^T^Pg, 
Pg = 1/ and for the diquark condensate with paired blue-green and green-blue ud 
quarks one has i] = 1/ y^. The penetration depth of the magnetic and gluomagnetic 
fields Xq and the mixing angle a are given by 

= ^^V^' cos a = . (4) 

At neutron star densities gluons are strongly coupled {q'^/Att ~ 1) whereas photons 
are weakly coupled (e^/47r = 1/137), so that a ^ r^e/gf is small. For g'^/4:TT ~ 1 we 
get a ~ 1/20. The new charge q is by an order of magnitude larger than e/Vs. 

Please notice also that since red quarks are normal in the 2SC and 3SC phases, 
there exist the corresponding normal currents = — II^'^A'^ and f^iG^j = 

—li^^^Ga which however do not contribute in the static limit under consideration 



to the above Ginzburg - Landau equations, cf . [T^ . Thus, the qualitative behavior 



of the static magnetic field for all three 2SC, 3SC, and CFL phases is the same. 
Recently the influence of a constant uniform chromomagnetic field on the formation 
of color superconductivity and the role of the Meissner effect for tightly bound states 



have been considered in the papers [119], pO|. 
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We will solve the Ginzburg-Landau equations (|l]), (||), (0) by properly taking into 
account the boundary conditions, in particular, the gluon confinment condition. We 
will find the distribution of the magnetic field in both the quark and hadronic phases 
of the neutron star and will show that the magnetic field penetrates into the quark 
core in the form of quark vortices due to the presence of Meissner currents. 

We assume a sharp boundary between the quark and hadron matter since the 
diffusion boundary layer is thin, of the order of the size of the confinement radius 



(~ 0.2 0.4 fm), and we suppose that the coherence length = y7/(— 2a) is not 
less than this value and the magnetic and gluomagnetic field penetration depth Xg 
is somewhat larger than the confinement radius. Also we assume that the size of 
the quark region is much larger than all mentioned lengths. 




2 Solution of Ginzburg-Landau equations 




superconducting matter 



(5) 
(6) 



where 



/ = Airiq-f [dVd* - d*Vd 



(7) 




(8) 



(9) 



then Eqs. and (|^) can be rewritten in the form: 




(10) 





rot TotA' + sin^ a A' 




sm a cos a 



We derive from equations (0) and (|T2D the relation 

rot rotGg = — cot a rot rot A' . 



(13) 
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We can consider instead of the system of coupled equations (1^) and (|TT]) the equiv- 
alent system (|T2|) and (|13|). We substitute Gg from (|12[) into (|1^) and obtain the 
following equation 

rot rot (A^ rot lotA' + A') = 0. (14) 



We introduce the new function M' as 



M' = rot TotA' 



(15) 



we obtain 



rot rotM' + M' 



0. 



(16) 



So we can determine the function A' by simultaneous solution of the equations ([T5|) 
and (0). Then we can find the electromagnetic potential A and the gluonic potential 
Gs from equations (|T2]) and (|^). 

In order to determine the distribution of electromagnetic and gluonic potentials 
inside the superconducting quark matter core in response to an external magnetic 
field we shall require on the quark-hadronic matter boundary both the continuity of 
the magnetic field and the vanishing of the gluon potential [Gg = 0) due to gluon 
confinement. As we shall see below, these conditions are sufficient for a unique de- 
termination of the magnetic and gluomagnetic fields inside the quark matter region. 

We shall assume that a neutron star with radius R possesses a spherical core of 
radius a consisting of color superconducting quark matter surrounded by a spherical 
shell of hadronic matter with thickness R — a. The functions M', A and Gs in spher- 
ical coordinates {r,{},ip) have only (/^-components M^(r, -i?), A^{r,{}) and Gg ¥'(^)f^)- 
For the solution of the equation (p6|) 



-ipV 1 " J 1 ^ '-'fix 

we make the ansatz M',„{r,d) 



Then equation ([16|) can be written as 



M^(r) 



sint?. 



r 



dr 



. 



(17) 



The solution of equation ([T7|) is 



M,(r) = - 



c; 1 



T \ J— 



Xn 



e^^ +Cn 1 + 



Air 



The condition that M^{r) tends to zero at the center of the quark core gives c'l = — c'a 
, so that 



M,(r) = ^ 



ly ly ly 

sinh -— cosh -— 

Xq Xq Xq 



(19) 



Substituting the solution ([T9|) into equation (|T^) for A' we obtain the following 
solution 

A'^{r,^) = M'^{r,'&)+CQ r sint? . (20) 
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Taking into account (^OD and using equation (||) for the electromagnetic potential 
we obtain 

A^{r,d) = M;(r,?9) +c[, r sint9 + . (21) 

The unknown function f^{r,i)) we will find from the solution of equation which 
gives 

/^(r,^) = Co r sint? . (22) 
Finally we get for the electromagnetic potential 



A^{r,i^) = M^(r) +Cor + 



Cor 



sin-i? . 



2 sin a. 

Using equations (|T^ and (^, we find the gluonic potential Gs^^ in the form 

Co r 



Gsifiir, {}) = — cot a M^(r) + tan a Cq r 



2 sin a 



sin?? . 



(23) 



(24) 



The constant Cq we can define using the gluon confinement condition on the surface 
of the quark matter core G8y,{a,{}) = . This condition will define Cq as 



, 2 ^y(Q) , Co 

Cq = cot a — 1 

a 2 sma 



(25) 



to be substituted into equations ( p3D and (p^ so that we obtain for the final expres- 
sions for the electromagnetic and gluonic potentials 



MJr) + cot^ a-MJa) + 

a sin a. 



sin-f? , 



M^(r) 



r 



-MJa] 



cot a sin '(9 . 



(26) 
(27) 



In ending this section we mention that the electromagnetic potential in the hadronic 
matter part of the neutron star can be found from the solution ([T8|) by replacing 
the penetration depth for quark matter Ag with that for hadronic matter Ap . 



3 The magnetic field components for the neutron 
star 

The components of the magnetic fields in quark and hadronic matter can be found 
from those of the vector potentials using the formula B = rot A . In spherical 
coordinates we have 



B^ 



d 



sint? A<^(r, 



r sin 'd di!) 
1 d 

— — (r A^{r,i^)) . 
r or 



(28) 
(29) 



For the case of quark matter we have to insert into these formulae the expression 
(refml) for A^{r,-^). Then finally we get (for r < a) 



r 



2M^(r) 



2 cot a- 



Mja) 



a 



+ 



-4- (r MJr)) + 2cot^a 
r dr 



2co ' 
sin a 

a 



cosi) , 

^ 2co 
sin a 



(30) 
(31) 



where M^(r 



is defined by equation (p^Qf). 
The magnetic field in hadronic matter phase we can find from the solution (|T8|) 
by taking into account that proton vortices in this phase generate a homogeneous 
mean magnetic field with amplitude B and direction parallel to the axis of rotation 
of the star For the components of the magnetic field B^ in the hadronic phase 
(for a < r < R) we get the following expressions 



Bl 



2AJr) 



+ B 



cos^ 



-^{r A^{r)) + B 
r dr 



(32) 

(33) 



where 

As we have mentioned above Ap is the penetration depth in hadronic matter. 

The external magnetic field B^ in the region outside of the star (r > R) has to 
be dipolar and their components are given by the following expressions 



Bt 



2M 



cost? , B' 



—rr sm^y 



(35) 



where M. is the full magnetic moment of the star. The unknown constants Cq, Ci, C2, 
C3 and M. in equations (p8|)-( ^ ) have to be defined from the continuity conditions 
of the magnetic field components at r = a and r = R and from the condition 



B''Vi + BV2 = 8n/3M 



(36) 



where -R'^ is the z component of the magnetic field in the quark matter region with 
the volume Vi. V2 is the volume of the hadronic matter region. Here we suppose 
that the magnetic field in both regions is mainly constant and parallel to the axis 
of rotation z. As we will see later this supposition is fulfilled with high accuracy 
because a, R and R — a are much greater than Ap and Xq. 

The continuity of the magnetic field components at r = a and r = R gives us 
the following equations 



2M^(a) 



2cot^^ 



+ 



2cn 



sma 



2AJa) 



a 



B 



(37) 
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1 d 
r dr 



+ 2 cot^ 



Mja) 



+ 



2co 



a sma 



R 



1 d 

r dr 



+ B 
+ B 



r=R 



1 d 
r dr 
2M 

'W ' 
_M 



r A^{r)) 



B 



(38) 
(39) 
(40) 



Substituting the functions M^(r) and A^{r) from equations (jTUp and ( PTf ) into the 
system of equations (|37|) - (^0]) we find the following system of equations 



Cl 



a \ a a a 
sinh---cosh- 



- C2 



1 _|_ '2 

Xp 3Ap 



/ a a 1 ^ 
Cl sinh — cosh — 



(41) 



sin^ a 



M 
M 



C2 [I 
C2 (1 



Xp SXj^ 
R\ JL 



R R' 



e^f +C3 I 1 + — + 
V 



R 

Arj 



An 3A^ 



A. 



e^p + C3 1 + 



R\ .JL 



A. 



\p _|_ 



p/ 
B 



where 



sin^ a — CqQ^ sin a . 



(42) 
(43) 
(44) 

(45) 



The solution of this system of equations using the fact that a, R and R — a are much 
larger than Xg and Ap we obtain the following expressions for ci, C2, C3 and 



Cl 



1 



aAp sin^ a + ^ sinh f- 



A. 



D e" 



C2 = 77^ 



2a sin^ a + 4^ sinh 



At, 



D e^p 



C3 



A„ 



2a sin2a + ^ sinh ^ 



B R^ R 



D 



1 



a sin" a + ^ sinh ^ 



(46) 

(47) 

(48) 
(49) 



In order to obtain a final result we have to determine the constant cq assuming B 
is constant. The latter is given since in the hadronic phase entrainment currents 



generate proton vortices which provide this field |21 
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Before we go over to the determination of Cq let us consider the behaviour of the 
magnetic field at distances r much larger than Ap and Xq. Also we take into account 
that A <C Ap -C a, A -C Ap -C i? and A <^ Ap -C i? — a. Using the expressions 
( ^BD - ( P^D for the constants, which we have obtained within this approximation, the 
components of the magnetic field in the different regions of the neutron star take 
the following form. 
For r < a 



for a < r < R 



2D 



a — r 



A„ 



sin^ a 



-u A 

"I" IT P 

At, 



2 cot a- 



2c, 







D 



( 

e ^1 



a — r 
An 



arAp sin^ a + 



+ 2cot^a- 



+ 

a sm a 



+ 



Co 



cosi) 



sin-i? 



a 



sm a 



(50) 
(51) 



and for r > i? 



BP 

r 



2D 



r — a 



ar^ sin^ a + ^ K 



B 



D 



r — a 



arXp sin^ a + 



B 



cos-i? , 



B 



cos-i? 



BR' . , 



(52) 
(53) 

(54) 



As can be seen from the obtained solutions, the magnetic field in both quark and 
hadronic phases depends on r only close to the phase boundary at r = a . In more 
detail, the part of the solution for the field which depends on r tends to zero in the 
quark core at distances a — r 3> Ag and in the hadronic shell at distances r — a ^ \p. 
But since Ap and Xq are much smaller than a and R — a except a small layer at the 
surface of the quark core with the depth of the order Ap + Ag, the magnetic field 
in both phases is constant and directed parallel to the rotation axis of the star, see 
solution (|50D - (^31). We conclude that in the main part of the volume of the quark 
and hadron phases the magnetic field is constant and has the direction z. In this 
approximation the condition ( ^6]) is satisfied. Therefore inserting in (p6D the relation 
M. = B R? /2, see equation (|^), we have 



5" 



2 cot a- 



Mja) 



+ 



2 Co 



B . 



Solving equations 



a sm a 
and ( |55D we finally obtain 
B sin a 



Co 



D = . 



(55) 



(56) 
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Thus in this approximation the magnetic field B enters from the hadronic phase 
into the quark phase in the form of quark vortices |TB[ . The transition zone is of the 
order Ap + \q which entails that the quantity D is small of the order (Ap + Xq)/a so 
that the condition D = is well fulfilled. 



4 Conclusion 

We have investigated the behaviour of the magnetic field of a neutron star with 
superconducting quark matter core in the framework of the Ginzburg-Landau the- 
ory whereby the simultaneous coupling of the diquark condensate field to the usual 
magnetic and to the gluomagnetic gauge fields has been taken into account (rotated 
electromagnetism). In solving the Ginzburg-Landau equations for this problem, we 
have respected the boundary conditions properly, in particular, the gluon confine- 
ment condition. We have found the distribution of the magnetic field in both the 
quark and hadronic phases of the neutron star and have shown that the magnetic 
field penetrates into the quark core in the form of quark vortices due to the presence 
of Meissner currents. 



Note added: 



After acceptance of this paper for publication in Astrofizika, the paper |22| has 
appeared on the e-print Archive, which complements the present study and adds a 
detailed discussion of the effects of rotation on vortices and gluon-photon mixing. 
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